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Abstract 

We consider a version of random motion of hard core particles on the semi-lattice 1, 2, 3, . . ., where 
in each time instant one of three possible events occurs, viz., (a) a randomly chosen particle hops to a 
free neighboring site, (b) a particle is created at the origin (namely, at site 1) provided that site 1 is 
free and (c) a particle is eliminated at the origin (provided that the site 1 is occupied). Relations to 
the BKP equation are explained. Namely, the tau functions of two different BKP hierarchies provide 
generating functions respectively (I) for transition weights between different particle configurations and 
(II) for an important object: a normalization function which plays the role of the statistical sum for 
our non-equilibrium system. As an example we study a model where the hopping rate depends on two 
parameters (r and (3). For time T — > oo we obtain the asymptotic configuration of particles obtained 
from the initial empty state (the state without particles) and find an analog of the first order transition at 
/? = !. 



1 Introduction 

In the famous paper [1] M.Fisher introduced models of one-dimensional random walk of hard 
core particles on the lattice. We shall consider a specific version of the models that Fisher 
called random turn walk models. They describe a motion of particles where at each tick of the 
clock a randomly chosen walker takes a random step. In both types of models each site may be 
occupied by only one walker at the same time. 

In our case we consider a version of this model where particles move along a semi-line, and 
where also a particle may be created at the origin. 

The model. Consider a set of nodes labelled by positive integers, where each nearest neigh- 
bors are linked by a pair of opposite arrows. Let us view the nodes 2, 3, 4, . . . as situated to 
the right of the origin, the node 1. An arrow which starts at a node i and ends at a node j 
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Figure 1: Configuration for T = 1 in case A*^^"°^ is the empty configuration. Black ball serves 
for a particle, white balls for free sites 



(where j = z ± 1) is assigned a weight equal to e~^^~^^\ Here U = (Ui, U2, . . . ) is a set of real 
numbers. 

As an initial state of our dynamical system related to a time T = 0, we place a certain 
number of hard core particles at nodes ("hard-core" means that each node may be occupied by 
at most one particle), this initial configuration will be denoted by X^^\ 

Consider a random motion of the hard-core particles where in each time instant one of three 
possible events occurs: (a) a randomly chosen particle hops along any of arrows attached to the 
corresponding node (i.e. either to the left or to the right) provided that target node is free of 
particles; (b) A particle is created at the origin (namely, at site 1) provided this site is free; 
(c) a particle is eliminated at the origin (provided that site 1 is occupied). We shall refer to 
configurations as neighboring ones if they differ by an elementary event from the list above. 

For instance, the empty configuration is an neighboring one to the configuration shown on 
figure [H because the last is obtained by the event (b) from the configuration without particles. 

The probability of each elementary (i.e. which occurs in one time instant) event is propor- 
tional to a weight (or, a rate) of the event which we assign as follows: 

(a) For the hop of a particle along an arrow i ^ j (where j = i ± 1) the weight is given by 
^-Uj+u^_ jj^g weight r(j) := e~^J+^% j = 2, 3, . . . will be called (right) hopping rate; 

(b) For the birth process the weight is (this weight will be called birth rate r(l)); 

(c) For the elimination process the weight is :^e^^ 

Along the random process in each each time step, say, T — > T + 1 a configuration A*^^) goes 
to a neighboring configuration A'-^^^-' with the following probability: 

where W^iT)^^ is the weight of the elementary event which creates (in one time step) a config- 
uration /i from the configuration \^'^\ the sum in the denominator ranges over all neighboring 
configurations. For instance, in case the initial configuration is the empty configuration, then 
by rules (a)-(c) with unit probability we obtain the configuration depicted in figure [TJ 

Each set of configurations \^^\ \^^\ . . . , A*^^^ where each pair A'^*\ A*^*+^^ is a pair of neigh- 
boring configurations will be referred to as a path of duration T which starts from the configu- 
ration A*^°^ and ends at the configuration A*^^^. The weight of the path is defined as the product of 
the weights of all elementary events along the path, Wx(o)^\wWx{i)-,\{2) ■ ■ ■ Wx(t--i)^xm- The 
transition weight between configurations A*^°-* and A'-^-' is defined as the sum of weights of all 
paths of duration T starting from A*^°) ending at A*^^^ and will be denoted by W^w^xi-t) (t). 



2 



Then, it is easy to see that the probability to come from a configuration A'^^-' to a given 
configuration A^"^) in T steps is the ratio 



W^A(0)^Aw(t) 
^A«.)(T) 



(2) 



where the denominator is a normalization function 



^A(0)(T)=$]W-A(0)^;.(T) 



(3) 



In what follows we shall omit superscripts for configurations. We shall not need intermedi- 
ate configurations which were introduced for better explanation of our model. 

Our letter is arranged as follows. In an introductory part to section |2] we introduce our 
tools: neutral fermions and the related Fock space, and quadratic operators (fTT]) - (fT2)) which 
depends on a given set U = {Ui, U2, ■ ■ ■)■ In subsection [XT] for arbitrary set U we shall obtain 
an explicit expression for the probability to achieve a given configuration in T steps in case 
the initial configuration is the empty state (the state without particles), see formulae (|27]) . (|34|) 
and (l35l) . For an arbitrary chosen external potential U it is impossible to obtain an asymptotic 
limit in formulae (|34l) and (|35]) in the large time limit. In subsection (|2.2I) we specify rate 
r(n) n = 1,2, .. . by formula (1381) . now the rate depends on two parameters: a constant 
denoted by r and an exponent p. In this case we present the asymptotic formulae for the 
density of particles, cr(n), see (|42|) and (|44l) . We shall show that for a given and large enough T 
all characteristics of the asymptotic configuration of particles: its size, the number of involved 
particles, the center mass e t.c. undergo a jump at /? = 1 (when T 00). We shall show that the 
normalization function Z)^(o)^q(t) = Z{r, /3, t) also has a jump. In our problem this function 
plays a role quite similar to the role of the partition function in statistical physics; in this sense 
in our model we may treat this jump at /3 = 1 as a first order phase transition. The expression 
for probability in different regions of the parameter /5 is given by (|55l) and (l56l) . Let us note 
that if we fix /5 = 1, then, our model turns out to be a discrete time version of the model called 
asymmetric simple exclusion process (ASEP), now the constant r plays a role of asymmetry 
parameter. In two parts of subsection 12.31 we shall link respectively transition weights (in item 
I) and the normalization function (see item 11) with tau functions of two different BKP, in the 
item II the key role is the relation (|64l) which relates rates r(n) to the BKP higher times. 

Let us notice that a wide usage of the free fermion approach to random partitions and certain 
random processes was presented by Andrei Okounkov in a series of papers, in particular see [2] . 
Our approach is different and based on (fT4l) . It is also different from the approach invented by 
H.Spohn and K-H.Gwa in [3] for the study of (the continues time) ASEP where spin system 
was used for the coding of particle configurations. In this approach the state spin up codes 
the filled stated, spin down the empty state. Via Jordan- Wigner tranform this spin system may 
be related to fermionic one. However in that case a quadric fermionic Hamiltoinian is used to 
describe the stochastic dynamics of particles which was identified with Hamiltonian dynamics 
of quantum spin (or, of nonlinear fermionic) system where the (real) wave function yields the 
probability distribution for the stochastic process. Following [4] we use free fermions and the 
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answer for the probability is given by the ratio of two factors (|35l) quite similar to what we have 
in thermodynamics where the probability of a state is the ratio of the weight of the state and a 
normalization function (the partition function). 

The main part of this work was reported by one of the authors (J.W.L) on the workshop 
"Random and integrable models in mathematics and physics" in Brussel, September 11-15, 
2007. 



2 Fock vectors and configurations of hard core particles 

Neutral fermions. In what follows we shall need neutral fermions, {(f)n}nez, as introduced by 
E. Date, M. Jimbo, M. Kashiwara and T. Miwa in [5], defined by the following property: 



[-iTK-m, (4) 

es me anucommuiaior. in particular, {(po''^ 

vacuum states by 



where [ , ]+ denotes the anticommutator. In particular, (0o)^ = ^- Next we define an action on 



0„|O) = O, (O|0-n = O, n<0, 0o|O) = -^|O), (O|0o = ^(O|, (5) 

Note that the definition of these "Fock spaces" is different from the usual one which was intro- 
duced in [5]. We follow [6], where the action of (po (which is not a number) is different. See 
Appendix for some more details. 

The basis of the corresponding right and left Fock spaces are formed by vectors 

|A) :=0,^...0;,^|O) (6) 

(A|:= (-1)1^1 (O|0_A^---0_Ai, (7) 
Ai > ■ ■ ■ > Aat > (8) 



and by 



where 



and |A| = Ai + \- Xn- 

We have one to one correspondence between configurations of hard core particles on the 
lattice 1, 2, 3, . . . and the basis Fock vectors. These configurations are also called Maya 
diagrams. Namely, the Maya diagram of the vector |A) is a set of vertices 1, 2, 3, . . . , where 
each vertex numbered by Aj is drawn as the black ball (a hard core particle), all other vertices 
are white balls. 

As we see, by dH)-© we have 

(A|/i) = 5a,m (9) 

Consider the following operator 

B = B^{U) + B^i{U) (10) 
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where U = {Ui,U2, ■ ■ ■ ), is a semi-infinite set of numbers (which may be also considered as 
variables), and where 

oo 
j>0 

(11) 

and 

oo 
j>0 

(12) 

It is straightforward to check the relations 

B_,B, - B,B., = ^ , 5_i|0) = = {0\B, , (13) 

which we shall need soon for the calculation of (|2T1) . At last let us note that for different purpose 
the operators Bi and were used in [7] (namely, to construct examples of multivariable 
hypergeometric functions which are also multisoliton BKP tau functions [8]). 

2.1 Stochastic system: description via fermions 

The birth-death of particles at the origin and their diffusive motion may be described as follows. 
A sequence of Fock vectors 

|A') ^ (B.iU) + i?_i(f/)) |A') ^ ■ ■ ■ ^ (B.iU) + B_^iU)y \\') ^ ■ ■ ■ (14) 

describes an evolution of the initial (basis) Fock vector | A') - where the variable T = 0, 1, 2, . . . 
plays a role of discrete time - to linear combinations of different basis Fock vectors. Due to 
the correspondence between basis Fock vectors and configurations of hard core particles, this 
evolution may be interpreted as the random process described in section [H where each time 
step is numbered by T. The details will be presented below. 
Let us notice that each basis Fock vector 

|A):=0Ai---0A..|O) (15) 

where A is a strict partition which is A = (Ai, A2, . . . , Aat), that is in one-to-one correspondence 
with a configuration of particles, located in the nodes with numbers Ai > • ■ ■ > Aat > 0. This 
is the reason, why we have chosen ^ as definition of the Fock space and not the one of [5]. 
For the Fock space introduced in that article this is not the case. 

We are interested in the discrete-time version of this random process which is given by (fT4)) . 
where each time step is numbered by T. 

As we see via ^ and Q, the first term in the right hand side of (fTT)) describes the creation 
of a hard-core particle located at the node number 1 (the origin) (provided that this node is 
free). Since 0o|O) = 75 10) we assign the weight -j=e~^^ to this creation process. Each other 
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terai of Bi describes a hop to the right. Similarly, the first term of i? i in the right hand side 
of (fT2l) describes the elimination of a particle located at the node 1 (provided that this node is 
occupied). Since again 0o|O) = -^lO) we assign the weight -^e^'^ to each elimination process. 
Other terms of i describes the hops to the left (backward motion in the direction of the 
origin). 

It is not difficult to see that 

Wy^xiU; T) := (A| {B,{U) + B_,{U)Y \X') (16) 

is a sum of weights of paths over all paths of duration T starting at the configuration A' and 
ending at the configuration A, and, therefore, yields the transition weight of the T-step random 
process which describes a transition from an initial configuration of the hard-core particles 
described by coordinates A'^^, . . . , A'^y, to a target configuration Ai, . . . , Aat defined in the In- 
troduction. As the number of particles does not have to be conserved along the process, N 
is not necessarily equal to A^'. Let us mark that for each path from a configuration A' to a 
configuration A of a duration T we have 

N -N' = n+-n_ (17) 

T = j++j-+n+ + n_ (18) 

|A| - |A'| = j+ -j_ +n+ -n_ (19) 

where j+ is the number of hops to the left during the time interval T, j is the number of hops 
to the right, n+ is the number of creations of a particle and r2„ is the number of eliminations of 
particles at the node 1. Let us notice that from (fTSi) and (fT9l) it follows that T and |A| — |A'| have 
the same parity. 

Consider the case A' = and look for the weight of the process which transports the initial 
empty state (there are no particles at all) to a given configuration A in T steps 

Wo^xiU; T) := (A| (5i(f/) + 5_i(f/))" |0) (20) 

In order to evaluate the right hand side we need the following formulae 

_ e"^-i|0) = |0) , (21) 

(the first equation follows from the Baker-Campbell-Hausdorff formula and from (fT3l)). and 
also the formula 

(A|e^^^(^)|0) = 2--.l^le-S£^^^. Hl^n Vfv ' (22) 

see [7], [8]. Here |A| := Ai + ■ ■ • + Aat is the weight of A. In our case |A|/A^ is interpreted as 
the location of the center of mass of the hard core particles. 

Now using (|2TI) . we obtain 

(A|e^^^(^)+^^-(^)|0) = e^(A|e^^^(^)e^^-(^)|0) = e'^(A|e^^i(^)|0) (23) 
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If we develop the left-hand side and the right-hand side of (1231) in powers of the variable z, we 
obtain the following formula for the transition weight from the vacuum state to a state A in a 
time duration T: 

T! / i \ ^_E. _Y"iv T— r i T— r Aj — Aj 



(A|(i?i(f/) + i?_i(f/))^|o) = iyo-.A(T) = ^ ( - j 2-f e-s-^^« nxrll 



Ai! . , . Aj + A,- 

(24) 



where the right hand side is non- vanishing only if the relation 

T=|A|+2m, m = 0,l,2,... (25) 
is valid. By (fTSi) and (fT9l) we obtain the meaning of m, viz.: 

m = j_ + ?7,_ (26) 

Thus 

[0 iffT-|A| odd, 

"'°-^<^'"> = |(^2l^l-fe-i:£.''^.nt.inf<,fe^ itfT-|A| even (27) 

The normalization function, counting weights for all possible target configurations, which 
may be achieved in the time duration T, is 

oo 

Z(t/;T):=^ E (A|(5-i + i?i)^|0) (28) 

iV=0 Ai>-->A]v>0 



°° iv ^ 2l^l-^ ^ e-^^> A A, - A,- 



(29) 



T— I A I eve.n 



Note that due to the Gamma function this sum is finite. 

The term corresponding to = gives in the summation the term (0|(i?_i + -Bi)^|0). 

Let us note that 

(1) |A| = corresponds to the returning to the initial position. This only occurs when the 
number of time instants is even, T = 2m, 

l^o^o(T) = y^2-^ = 2-5(T- 1)!! (30) 

(2) Now consider the next case where the final configuration is the one-particle one. Now 
A is a number which denote the coordinate of this particle. For simplicity we consider rin) = 
r = const. Then, given large enough T, by Stirling's formula we can present the transition 
weight as follows 
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where 



Wo^3^ = ■ ^r^-'e-""', Ai = rV2T-r^-- + 0(T 

r(A + i)r + 



The formula resembles formula for the Brownian motion (here the variance is given by T4). 
At last we note, that in case the rate r depends on site, then, in a wide class of rates in large T 
limit Ai may be evaluated as the solution of Ai = r{Xl)^/2T. For instance, for Gauss potential, 
Un = \n'^, one obtains Ai ~ log(2T). 

(3) For given A in large T limit (this means that T ^ |A|) by the Stirling's approximation 
we have 

Tyo^^(T) = H^o-o(T)T^2-^e-^2l^le°(^°)e-^\ (31) 



where 

N 

^--lo^nvn^ (32) 



. -, Aj! . Aj + Aj 

means an electrostatic energy of Coulomb particles (placed in an external field) which are 
attracted by their image. We see that in the large time limit the weight of a configuration 
increases with |A|, and for given T and | A| depends only on the Ex. 

(4) The case T = |A| = n+ + j+ corresponds to the non-stop creation + forward motion 
processes. Then 

W,^,{T) = T!2-f e- ^^^^ n ^ n (33) 

Aj. Aj + Ao 

1=1 t<j ■' 

In case the potential is a rapidly decreasing functions f/„_i » Un (and, therefore, left hopping 
rates are much larger than right hopping rates), then the configurations where T = |A| are 
dominant in the sum for normalization function. 

Let us note that up without the factor e~ ^'=1 the number I^o^a(t) is equal to the 
number of shifted standard tableau of shape A, see [16], that is the number of ways the Young 
diagram of the strict partition A may be created by adding box by box to the empty partition in 
a way that on each step we have the diagram of a strict partition. 

Finally, one can get rid of the restriction Ai > ■ ■ ■ > A at in the summation (|29l) rewriting it 
as a sum over all non-negative integers Ai, ■ ■ ■ , Aat: 



Ar=0 Ai,...,Aiv>o -^12 ^ ^) i=l * i<j 



T— I A I even 



Aj + Aj 



(34) 



The probability to come to a configuration A in T steps starting from the vacuum one is 
given by 

-^(--'-^ 

In the present paper we do not write down the expression for the probability p^^xiU] t) 
because it is rather cumbersome and contains the so-called skew projective Schur functions. 
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2.2 Asymptotic configuration of the particles in T ^ oo limit 

One may ask what configuration is obtained from the empty configuration (vacuum) configu- 
ration due to the creation/annihlation processes at the edge of the lattice and to the hops of the 
particles in the large T limit. (Such configuration will be called asymptotic one and denote as 
A(t)). To find it we should find the largest term in (|34|) . Let us do it via saddle point method. 
First, in the large T limit it is reasonable to introduce the density of particles a{u), where u = ^ 
where R is the size of configuration. The number of particles is 

N = R [ a{u)du (36) 
Jo 

Let us note that the density interpolates between full package state ((t(m) = 1), and empty 
state {a{u) = 0): 

< a{u) < 1 (37) 

Let us find the asymptotic configuration in the case when the creation rate is r (1) := = 
r and the external potential is 

Un = — nlogr + {f3 — 1) \ogn\ , n = 2, 3, . . . 

which means that creation and (the right) hopping rates are as follows 

r{n)=rn^-^, n= 1,2,3,... (38) 

P > 1 describes a locking potential while in /5 < 1 case the potential try to drive particles to 
the right from the origin. The case (3 = 1 may be considered as a discrete time version of the 
so-called asymmetric simple exclusion process (ASEP) on the half-line, now, the parameter r 
being an asymmetry parameter. 

1 

Let notice that formally the point n = = r^-i^ h di point of an extremum of the potential 
Un where the left and the right hopping rates are equal. 

It means we want to find the configuration A = A(t) where for given T the weight is as 
follows 

Wo^x{r,(3-T) = T!-^^2-f e-^^^^''^) (39) 
here T — |A| is even, and the"electrostatic energy" of the configuration is 

S.(r./J)^-.og(rNn^n|^) (40) 

In the continues limit in a standard way one can write the saddle point equation for sum 
Z{t) which will define the density function a in the large time limit, T ^ oo. For A G (0, R) 
we get 

(41) 
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where P f stands for the principal value. This equation is to be solved by a standard method 
[22], see Appendix. For < /3 < 2 we find 

a{XR-^) = - arccos ^, A G [0, R] (42) 

TT it 

(The validity of (l42l) may be verified by substitution in iAll). The constraint (1371) causes the 
restriction < /? < 2 for our solution in formula (|42)) . 
The weight of such configuration is 

|A(T)| = f ua{u)du = ^ (43) 



and after substitution into the logarithmic term in (|4TI) we obtain a relation of i? to T as follows 

3 22"^^ 



Let us make a remark. In the approximation we consider one can replace sums by integrals only 
up to terms of the order 0{R^). It results to a fact that r in the relation (|44l) should be replaced 
by an effective (yet undefined) hopping rate r^^f- = t/tq, where tq is of order e*^*^^"^. Below 
by r we will imply this effective hopping rate. Arguments that r^^ = r will be considered 
separately. 

Equation (|44)) shows that in the large T limit the dependence of i? on T is different in regions 
0<P<1, P = l and 1 < /3 < 2. Let us introduce 

R, = 2(/5r2)^, T, = /3{/3r^)^^ 



As we see in the region when R = R^, the terms in the r.h.s. of (I44|) are equal. Below we imply 
that T >> T*. In the large T Umit as we see from (|44|) 



R = R{p,T) = < 



8t 



1+r 



-2 



if < /3 < 1 

if /5 = 1 (45) 



2 {2r^j) 2« if 1 < /? < 2 



Thus, when < P < 1 the size of the asymptotic configuration is proportional to y/r, while 
in the vicinity of /3 = 2 we have the forth root behavior. As we see the discontinuity appears 
at /? = 1 in the large T limit. The same behavior has the number of particles in the asymptotic 
configuration which is proportional to the size of the configuration: 

N{r, (3,T)=R j^^ arccos du = ^^^lll (46) 

The weight of the asymptotic configuration (|43l) for large enough T is 

2 't + 0(t^) ifO</3<l 
|A(T)| = ^^<!^ ^ if (3=1 (47) 
f(2r2T)^ ifl</3<2 
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Figure 2: The evolution of the empty configuration in T ^ cxd limit for /5 = 1. Asymptotic 
density of particles. 



Notice that the center mass of the configuration also has a jump at (3 = 1: 



|A(T)| 
iV(r,AT) 



TT 



8(l+r-2) 



if < /3 < 1 
if (3 = 1 



(48) 



f {2r^T) if 1 < /3 < 2 



L 4 



For large enough T, the number m(T) = j-(T) + ^^-(t) depends on region as follows 

'O(t^) 



m(T) = -(t- IAI) = 1t-^ 
^ ^ 2^ ' 2 16 



if < /3 < 1 
if /3 = 1 



T_|(2r2T)^ ifl</5<2 



2(l+r2) 

T _ ^ 
2 4 



(49) 



In (|47l) and (1491) we keep terms which we shall use in evaluations below. 
For large enough T, the electrostatic energy (l40l) of A(t) is 



EA(T)(r,/5) = -|A(T)|lnr + 

pR^ [ (T{u){u\n(uR) -u)du^ --R^ [ [ (T{u)(T{u')ln^^^^dudu' (50) 
Jo J ^ Jo Jo u + u' 



|A(t)|(/5 \nR-\nr)- R'(3'A 



where 



j3 A = —(3 / lyiu) ivXmi — u) du + / du' / cr(u)cr(u')ln du = (3 

' ' ' u + u' 



u — u 



1 In 2 
16 ~8~ 



(51) 



see Appendix IA.4I which yields A 
By (05]) and (07]) we obtain 



_L _|_ In 2 
16 8 • 
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f TliiT - T (f Inf + f + In r) if < /3 < 1 

^A(T) (r, (3) ^ { + t:^ {>u^ - 8A - In r) if /5 = 1 (52) 

I (2r2T)^ (InT + (1 + 2(3)\n2 - 16/3A) if 1 < /3 < 2 

Using Stirling's approximation to evaluate m(T)!, for T large, we obtain the weight of the 
process as follows 

o— 2m(T) ^, , 

iyo^A(T)(r,/3;T) = t!^— 2-^e-^MT)M) (53) 
m(T)! 

'^,g-|TlnT+T(|ln| + |+lnr)+0{A) if < /5 < 1 

Ti„ I ■pr ln2(l + r2) , ln2(l + r-2) , 21n2-l + lnr2 , 1-ln 4 i , 
^lg-2i°T+T| + 2(^_^^_2) + 2{l+r-^) + 2(l+r^) >+'^^VTJ ^ = 1 (54) 



^iyo-o(T)e ^ ' y J ifl</5<2 

Then, as we see the normalization function Z{t) which according to the saddle point 
method has the same leading term in the large T limit as 11/o^a(t) (r, (3] t) has a discontuinity at 
(3 = 1 which may be interpreted as a sort of the first kind phase transition in our non-equilibrium 
system. 

Now we can evaluate the type of asymptotic of the probability to achieve a given configu- 
ration in T ^ oo steps. We have Z{t) = 11/o^A(T)(T)e*^''^°^\ where the last factor originates 
from the Gaussian integral around the saddle point configuration A(t). Then for T 3> |A| we 
have 

Po-.A(r,/3,T) ^ ^ T^e-^2l^le-^^('-'^)e-(^'^'^) (55) 

WQ_x{T){r,p,T) 



where u does not depend on A. (For the enumerator of (p5l) we used (1311) and Stirling's approx- 
imation for (l30l).) 

The answer depends on the region of (3 



e 



a;(T,r,/3) 



^/3^TlnT-T(|ln| + |+lnr)) + lln{2e)+0(v^) jf Q < /3 < 1 
}+0(v^) 



f (2r2T)'3(l-/3)+0(T2? 



2 "^2(n-r-2 
i / 1 



if /? = 1 (56) 



e"' ' V y if 1< /5 < 2 

where h = 21n 2 — 1 ^ 0.4. As we see, in each case, in the large T limit e"^ is vanishing. 

At last let us note the following. As we see in case of a decreasing potential (or, the same, 
a increasing rightward hopping rate), < /5 < 1, the weight of the asymptotic configuration 
is equal to T which means that the asymptotic configuration is created by only creating events 
at the origin and rightward hops, there were no elimination events and backward hops in the 
history of this configuration, j_ = n_ = 0. 

For /? < solution does not exists. In /3 — >^ +0 limit the number of particles (|46|) vanishes, 
while |A(t)| is equal to T. Indeed, the external potential f/„ is decreasing so rapidly that the 
largest weight has the one particle configuration where the particle moves in the ballistic way: 
it is located at the distance T to the origin. When /3 > 2 we have a locking potential which 
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forces particles to form a sort of a condensate near the origin where all sites are occupied. 
The size of the condensed phase is defined by p. This problem is treated by a method similar 
to the suggested in [21]. Along this way we can show that the solution is given in terms of 
elliptic integrals of the first and third kind. The problem will be considered in detail in the next 
version of this paper, where we will show that the normalization function Z(t) has singularity 
as function of the parameter P at P = 2. 

2.3 Relation to the BKP tau function 

The link of the described stochastic system to integrable equations is two-fold. 
(I) A BKP tau function as generating function for transition weights Wy^xir). 
Let us consider the following vacuum expectation value 

t{s,s;U,z) := (0|e^(^)e^^e^(^")|0) (57) 

where B is given by (flOl) and 

ri=l,3,5,... meZ n=l,3,5,... meZ 

(58) 

with z, U = {Ui, f/2, • • • ) and sets s = (si, S3, S5 . . . ) and s = (si, S3, S5, . . . ) parameters. 
Function r(s, s; U, z) depends on f/ = (f/i, f/2, . . . ) as the operator B depends on these param- 
eters. 

The hierarchy of of Kadomtsev-Petviashvili equations of type B (the BKP hierarchy) was 
introduced in [5]. As we have already mentioned we use its modification suggested in [6]. 
These is a semi-infinite set of compatible nonlinear differential equations which may be viewed 
as a set of commutative time flows. It is common to enumerate the BKP equations (flows) by 
odd numbers. 

Expression (l57l) where U may be chosen as an arbitrary set of numbers, provides an example 
of the BKP tau function constructed in [6], where s = (si, S3, S5 . . . ) is the set of the so- 
called higher times (the times of commutative flows related to different equations of the BKP 
hierarchy). Actually the set s = (si, S3, S5, . . . ) is also related to a (second) BKP hierarchy of 
equations, which is compatible with the first one; thus, (157]) is an example of the tau function 
of a coupled BKP hierarchy. 

Using results of [9] and doing similar calculations in the framework of BKP hierarchy con- 
structed in [6] (see Appendix) one can show that the operator e^*^*^ applied to the left vacuum 
generates all basis left Fock vectors, quite similar e^^"*^ applied to the right vacuum vector gen- 
erates all basis right Fock vectors. Namely, we have the following left and right coherent states 
(see also ©) and ©) 

00 

(0|e^W = $:2"f Yl ^^(0(^1' (59) 

N=Q Ai>-->A]v>0 
00 

e^(^)|0) = $:2-f Y: ^^©1^) (60) 

Ar=0 Ai>-->Ajv>0 
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where Qx (|) , A = (Ai, A2, . . . ), are known as projective Schur functions, see [16]d, which, in 
an appropriate space, form a complete set of weighted polynomials in the variables si, S3, . . . . 
Using these formulae we obtain that the BKP tau function (|57]) is the generating function for 
(fT6l) . namely 

00 00 -J- _ 

ris,s;U,z) = Y, E '^'^''^ E ^^3a (0 (|) W^v^A(f/; t) (61) 

T=0 N,N'=0 Ai> ..>Ajv>0 

What one obtains in case of a general BKP tau function will be explained in a more detailed 
forthcoming paper. 

(II) A tau function of dual BKP as the generating function for the normalization func- 
tion Z(t) 

Let us introduce the generating function for the normalization function (l28l) as follows 

Z(f/;z):=E E (A|e^^|0) = E^^(f^;T) (62) 

Ar=0 Ai>--->Ajv>0 T=0 

By (EU) and (ES) we have 

N=0 Ai>-->Ajv>0 i=l i<j * ^ 

--E^ E n^n 

N=0 Ai,-- ,Aiv>0 i=l i<j 

Let us mark that the last sum may be compared with the grand partition function of the log 
Coulomb gas, compare with [18], [19]. 

Now, we want to present parameters U = {Ui,U2, ■ ■ ■) as follows 

00 

Un = - E (64) 

m=l,3,5,.-- 

where t = (ti, ts, t^, . . .) is a new set of parameters (the same parametrization was used in [8]). 

Remark 2.1. Let us notice that our main example (|38] ) is not well fitted into such parametrization. It 
seems that one needs to introduce an additional flow parameter 

00 00 
Un = - ^ n"^tm + /?ln n! = - ^ n^tm + Pin In n-n) 

m=l,3,5,.-- m=l,3,5,... 

and consider the discrete dynamics with respect to /?. In the large T limit we may replace In n! by 
n In n — n, then it may be related to the flow introduced in in [14]. 

"^Here we use notations of [9]. See also [10] for the relation of the projective Schur functions to the BKP 
hierarchy. 



Ai — Ai 



Aj + Ai 



(63) 
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Then, the right hand side of (|63l) may be written as the following vacuum expectation value 
(compare with [11]) 



Z(t/;z) = e"r(0|e^We^(^) 1 + 72^]— 0H0o |0) 



(65) 



n>0 



.n+m 




(66) 



The proof is similar to the proof in [1 1] its sketch may be found in the Appendix lA. 51 

Now we notice that (|65l) is another example of the BKP tau function [6] which may be 
related to the so-called resonant multi-soliton solution where the number of solitons is infinite 
and momentum of solitons are nonnegative numbers (compare with [19]). Thus, the tau func- 
tion of this (dual) BKP hierarchy is a generating function for normalization functions. Via (|64|) . 
the higher times t = {ti,t^,t^, . . .) of this dual BKP hierarchy parametrize hopping rates of 
the particles of our stochastic model. 

Let us we note that time variables t = (ti, ^3, ts, . . . ) are integrals of motion for the BKP 
hierarchy mentioned in (I), in this sense the BKP (II) may be referred as a dual to the BKP (I). 

At last let us note that the example (l38l) where we put r = e*^ is related to 



One may conjecture that Z(ti,P;z) is a solution to a difference (with the respect to /?) - 
differential (with respect to ti) Hirota equation. 

Discussion 

In the continuation of this paper we shall consider our model where the parameter /3 > 2 
describes a model where a condensate of particles (the region of full package) fills a region 
near the origin. We shall explain the appearance of a phase transition at /3 = 2 which is 
rather similar to the transition studied in [21]. The other model where the injection rate is a 
free parameter will be considered where an analog of the phase transition will be presented, 
then it may be interesting to discuss links with [15]. It is interesting to understand links with 
results of [20] and with the approach of [18]. It may be also interesting to link our results 
with the results of a recent paper [24] where in the context of the quasiclassical limit of Toda 
lattice hierarchy of integrable equations the Vershik-Kerov limit shape for random partitions 
was reproduced. 




.n+m 
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A Appendices 

A.l A remarks on BKP hierarchies [6] and [5] and related vacuum ex- 
pectation values 

Let us note that different vacuum states were used in the constructions of BKP hierarchy in 
versions [6] and [5]. If we denote the left and right vacuum states used in [5] respectively by 
'(0| and |0)' then 

(O| = ^'(O|+'(O|0o, |0) 

Introduce also 

'(l| = v^'(O|0o, 
then'(0||0)' = = 1 and instead of © we have 

0„|O)' = 0„|1)' = O, '(O|0_„= '(1|0_„ = O, n<0 (69) 

see [5] for details. 

Correspondingly Fock spaces used [6] and [5] are different. From the representational 
point of view this definition is somewhat more convenient, since each Fock module remains 
irreducible for the algebra Boo which is the underlying algebra for KP equations of type B 
(BKP), see [6]. 

The vacuum states '(0| and |0)' are more familiar objects in physics. In particular any vac- 
uum expectation value of an odd number of fermions vanishes, while, for instance, (O|0o|O) = 
1 

Let F be a product of even number of fermions. Then it is easy to see that 

(0|F|0) = '(0|F|0)' (70) 

Let us note that all vacuum expectation values (v.e.v.) used in the present paper (say, tau 
functions (|57]) and (l65l) ) are sums of v.e.v. of monomials containing even number of fermions. 

In the context of applications to random process we can construct a model of random motion 
on the semi-infinite lattice based on the Fock space used in [5], however we find that the Fock 
space used in [6] is much more natural to our point of view. 



V2 

v^0o|O)' 



(67) 



(68) 
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A.2 A remark on formulae containing Qx functions 

Here we shall show that formulae found in [9] are of use in our case (see also [8]). 
From [9] it is known that 

for A/ odd, 



V2'(O|0oe^(^VAi0A,---0A^|O)' 



2-^Q(A„A„...,A^)(t) foriV odd, 
for even, 



where Q\ (|) , A = (Ai, A2, . . . ) the projective Schur functions, see [16]. Thus 

(O|e^(^)0A,0A,---0A^|O) = 



1, 



(O|e^(^VA,0A, ■ ■ -^A^IO)' + '(O|0oe^(^VA,0A, ■ ■ ■0A^0o|O)'+ 



1 

3 



+ ^ '(O|0oe^(^VA,0A. ■ ■ ■ 0A^|O)' + ^ '(O|e^(^)0A,0A, ■ ■ ■ 0A^0o|O)' 



1 

7i 



1 

2 



+ -^'(O|0oe^(^VA.</'A. ■ ■ ■</)A,|0)'+ i='(l|0oe^W0A.</>A. ■ ■ ■</>A.|l)' 

2"~Q(Ai,A2,...,Aiv) 

since the role of '(0| and '(1| (resp. |0)' and |1)') is interchangeable. 



(71) 



(72) 



A.3 Case < < 2: Vershik-Kerov type asymptotic 

We want to solve the following singular integral equation where r is a given constant: 



X — X 



P 



a{xR ^)dx 1 
X + x ^ 2 



ln2 T 



R 



xa{xR'')dx = (73) 



where X E [0, R] and where a is defined on the interval x E [0, R]. Here and below P serves to 
denote the principal value of integrals. 

Remark A.l. Let us notice that the second term in ([73] ) describes the repulsion of charges distributed 
with density a along {0,R), while the third term may be interpreted as the attraction of these charges to 
their image in the mirror. Then, it is natural to continues a to the interval {—R, 0) such that <j{—x) := 
— (j(x) (this describes the replacing of particles (with coordinate, say, x) by holes (situated in —x). 

p a{xR-^)dx _p a{xR-^)dx _ p a{xR-^)dx 
Jo A - X Jq X + x J_ji X-x 
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Notice that in this case in general we get a jump of o" in x = 0. For our future purpose (of in- 
verting Hilbert type transforms) we prefer to deal with continues functions. For this purpose we shall 
modify a in the region {—R, 0) by adding a constant equal to this jump, A, getting continues modified 
a (compensating this change of a by adding logarithmic term to the right hand side), see below 

Now, a{x) is a function defined on the whole interval [—R, R] via relation 

a{-x) = A - a{x), x E [-R, R] (74) 



(As we shall see later A = /3 if < /3 < 2, and A = 2 if /3 > 2). Then , we re-write U3 in a 
more convenient form as 

P r "^^""^^"'^"^^ = ln(/? + A)'^ + InA'^"^ - In(Cr), A G [-R, R] (75) 

J-R A — X 



-R 

where we denoted 



:= 2 ( T - / xa{xR-^)dx ] (76) 





C is an independent of A constant to be defined later. 

Singular integral equation (l75l) may be solved by standard method, see [22]. In case 
a{xR~^) is continues and bounded on [— i?, R\ the solution is given by formulae 

a{xR-') = -7r-2p / In^ -, x G [-R,R], (77) 

J~R V A2 - i?2 Cr X- X 

{5 := P — A) see formula (42.26) in [22]. Let us evaluate the integral in an explicit way. 
We consider y/X^ — R"^ in the integrand of (1771) as single- valued function with the cut [— -R, R] 
whose upper limit on [—R, R] is positive and lower limit on [—R, R] is negative. Also we 
shall consider ln(A + i?)"^ in the integrand as the single valued function defined on the whole 
complex plane with the cut [— oo, —R]. The cut of InA^ will be viewed on the ray {—oo, 0), a 
little bit above the real axe. Then we have 



Vx^-R^ (A + R)^X^ _ 1 r Vx^-R^ (A + R)^X^ dX 

J^R VA2 -R^"" CV l^~~2Jc7W^'' Ci- 

where the contour C is going crossing the points and 2R: C = — C_|. + C_. One can inflate the 
contour through the point 2R as there are no 'bad' singularities there. We have a cut (— oo, 0) 
caused by the logarithm. Inflating the contour to the infinity we see that the only contribution 
will be caused by the cut of the logarithm (we come to the contour which starts on minus infinity 
going a little bit upper the real line, than turning at the origin and going back to minus infinity 
a little bit under the real line. (The contribution of the circle embracing infinity vanishes as the 
integrand's asymptotic is |^.) This yields the integral along the cut 

^ V^^^. (A + i?)^A^ dX 



1 /-^ - i?2 dX 1 f° Jy^-2Ry, , dh 

' .{27TdV^) r--/ / ^ (27rAv^) (80) 



2 J-R VA2 - i?2 'x-X 2 ^h^ - 2Rh 'y-h 
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where h = \ + R, y = x + R and 27r A-\/^ and 2nA\^^ are the jumps of the logarithm. 
The first integral is [23] 

ZTTOV— 1 T = arcsm 



R 



VA2 - x-X 2x 



As we see taking small x this integral is not a real number; this results in the condition S = 0, 
or, the same 

A = (3 



For the second integral of (1801) we have (see [23]): 



^ ^ (y ^\ 

arccos 1^ — 1) (81) 



VWT2Rhy + h ^-y2 + 2i??/ vi? 
Inserting the last formula into (|80l) and then into ([771) we finally obtain 



a(xi?"^) = - arccos ^, x G [0, i?], < /? < 2 (82) 



We add the last un-equality to provide (1371) . < cr < 1. 

At last let us mark that one may show that if we relate the asymptotic configuration (|82|) 
to a strict partition, then, the shape of the Young diagram of the double of this partition will 
coincide with the so-called Vershik-Kerov asymptotic shape [12], [13]. 

A.4 Two useful integrals 

One can show that 



In 2 



1 /"^ 1 

— / arccos u {ulnu — u)du = — 

71" Jo 8 
^ '•^ u-u' , 1 



(83) 



— / du I arccos u arccos u' In du' = (84) 

TT^ Jo Jo u + u' 16 

which yields A = ^ + 

A.5 The relation between formula d65]) and (1631) 

Recal from (l66l) that (p^z) = ^^'Pi- Using dH) and ([5]), we can also calculate the vacuum 
expectation value for two of these fields: 



m>0 

(85) 

2 zi + Z2' 
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(where we assume that \zi\ > \z2\). Using Wick's Theorem we obtain that 



(O|0(Zl)0(;^2)■•■0(^2„^)|O) 



(O|0(zi)0(2;2)---0(^2™+i)0o|O) 



n 



l<i<j<2m 
m+1 



Zi Zj 
Zi -\- Zj 



(86) 



n 



Zi Zj 
Zi -\- Zj 



l<i<j<2m+l 

For the exponentials of the Hamiltonians H{t) and H{t) introduced in (l58l) one has 

e^W|0) = |0) (0|e^W = (0| 

and 



where 



hit, z)<i){z), 

h{t, z) = exp ^ tnz'^ J 

\n>l, odd / 

Using dHIl), dHT]) and ([88]) we calculate 

m 

2^ (O|e^W0(^i)0(a:2)---0(x2„^)|O) J]s^^n(a;2.-i-X2. 

1=1 



(87) 
(88) 
(89) 



all Xi>0 
2m 

allxi>Ofc=l l<i<j<2m * ^ i=l 

2m 

E E 11"^^.) n 

7rG52m a;^(2m)<a;,r{2m-l)<---<a;^(i) k=l 



Y\ sign{x2i-i - X2i) 



2m 



: 5^ 5^ stgn{7r-')l[b{t,Xk) J] 

Te52m 3;2m<a;2m-i<---<a;i 



k=l 



l<i<j<2m * -^1=1 



- n s^£/^(^~^(20 - 7r"^(2^ - 1)) 



2m 



l<«<i<2m * 



2™m 
(2m) 



3::2m<2:2m-l<---<2:i fc=l 
2m. 

f E n^(^^^) n 



all Xi>0 k=l 



l<i<j<2m 



H~ 



(90) 



We have used that 



Y sign{TT-^)Ylsign{Tr-\2i)-n-\2i-l))=2"'m\ (91) 

which is obvious if one realizes that the right-hand side of (|9TI) is equal to 2"^m! times the 
Pfaffian of the 2m x 2m matrix A = (Aj)i<j j<2m, where Aij = sign{j — i). 
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Next, we calculate 

2^+1 J2 (O|e^W0(xi)(/)(x2) • • • 0(x2^)0(x2^+i)0o|O) n st9n(x2i-i - x^i) 

all Xi>0 i=l 
2m+l m 

= ^ Yib{t:Xk) Yl ' [ Ylsign{x2i-i - X2i) 

all Xi>0 k=l l<i<j<2m+l ' ^ i=l 

2m+l m 

^ Yl n n ';;ri~;^]1^^9n{x2i-i - X2i) 

7re52m+l a;„(2m+l)<^7r(2m)<---<3;,r(l) ^=1 l<i<j<2m+l * 1=1 

2m+l m 

= E E sign{7r-')l[b{t,x,) J] J] ^^^^(^"'(2^ - ^"'(2^ - 1)) 

■"■£52™,+! 3;2m+i<a;2m,<---<a:i fe=l l<i<j<2m+l ' i=l 

2m+l 2m+l m 

= E E E ^^^^(^-') n ^(^' n n «^^^(^-'(2^) - vr-^(2z - 1)) 

(=1 TeS2m+l, X2m+l<X2m<---<Xl k=l l<i<j<2m+l * 1=1 

■7r(f)=2m+l 

2m+l m 

= E E ^^^^(z^) n ^(^' n n ^^^^(p(2o - p(2^ - 1)) 

/oeS2m a:2m+i<a;2m<---<a:i fe=l l<i<j<2m * i=l 



2m+l 



X2m+l<X2m<---<Xl k = l l<i<j<2m+l * 



E n^(^'^^) n 



(2m + 1)! 

^ ^ all Xi>0 k=l l<i<j<2m+l 



X I 3/ 



(92) 



Here we have used that 



Y sign{n~^) sign{7r~^{2i) - 7r"^(2i - 1)) = 



"€■32771+1' «=1 
7r(£)=2m+l 



-(-)' E sign{u)Y[sign{u{ee{2i)) - u{ei{2i - 1))) ^ 



"£■32771+1. i=l 



-{-)' sign{p)\{sign{p{2%)-p{2i-\)), 

P&S2m i=l 

where €e{i) — iifi < £ and €e{i) = i + 1 if i > ^. So if 

_i _ A 2 ■ ■ ■ £ - 1 i £ + 1 ■ • • 2m 2m + 1 

Vii h ■■■ je-i ji je+i ■■■ hm ^ 

then 

'1 2 ■■■ £-l££+l£ + 2--- 2m 2m + 1 

\3i 32 ■ ■■ je~i ^ je je+i ■ • ■ 32m-i 32m 
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and 

fl 2 ■■■ ^-l t £+1 ■■■ 2m -1 2m\ . fjfc if Jfc < 
p = . . . . . . . , «fc = < 

For instance for £ = 3 and m = 2 one has 



1234 _/1 234 5\_^ _A234 

5 1 4 2 3J^'^~l 5 1 3 4 2j~*^~U 1 3 2 



Hence, if we develop (1651) as a power series in z and keep the term e 4 . Then substitute (190] 
and dH) and use dMl), we obtain (l63l) . 
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